ACEEE Int. J. on Control System and Instrumentation, Vol. 03, No. 02, March 2012 



Optimal Feed-back Switching Control for the UPFC 

Based Damping Controllers 

Yathisha L 1 and S Patil Kulkarni 2 
'Research Scholar, E& C Dept, S J College of Engineering, Mysore, India. 
Email: yathisha_171 @yahoo.co.in 
2 Associate Professor, E & C Dept, S J College of Engineering, Mysore, India. 
Email: pk.sudarshan@gmail.com 



Abstract — This paper presents an optimal feed-back switching 
concept for the Unified Power Flow Controller (UPFC) based 
damping controllers for damping low frequency oscillations 
in a power system. Detailed investigations have been carried 
out considering switching between two optimal damping 
controllers; one with respect to modulating index of shunt 
inverter ni^ and another with respect to modulating index of 
series inverter nig- The proposed UPFC switching model 
presented here is tested on the modified SMIB linearised 
Phillips-Heffron model of a power system installed with UPFC 
using MATLAB/SEVIULINK® platform. The investigations 
reveal that the proposed optimal feed-back switching control 
between UPFC damping controllers nip and tllg provides 
moderately better performance with respect to settling time 
for both individual controllers as well as coordinated damping 
controller. 

Index Terms — OFSC, COC, UPFC, LQR, SMIB, Phillips- 
Heffron Model. 

I. INTRODUCTION 

The Unified Power Flow Controller (UPFC) is a multi- 
functional flexible AC Transmission (FACTS) device, whose 
primary duty is power flow control. The secondary functions 
of the UPFC can be voltage control, transient stability 
improvement, oscillations damping. It combines features of 
both Static Synchronous Compensator (STATCOM) and 
Static Synchronous Series Compensator (SSSC). 

Design of control strategies using FACTS devices such 
as UPFC for optimal power flow with improved performance 
is a major research concern of power system control 
community. Wang [1] has presented a modified linearised 
Phillips-Heffron model of a power system installed with UPFC 
and addressed basic issues pertaining to design of UPFC 
based power oscillation damping controller along with 
selection of input parameters of UPFC to be modulated in 
order to achieve desired damping. Wang has not presented a 
systematic approach for designing the damping controllers. 
Further, no effort seems to have been made to identify the 
most suitable UPFC control inputs, in order to arrive at a 
robust damping controller for optimal performance of all the 
state variables. However, in recent times, researchers are 
working on the selection of UPFC control parameter for the 
design of UPFC damping controller by applying different 
control techniques like Phase Compensation, Fuzzy Logic, 
optimal control techniques like Linear Quadratic Regulator 
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(LQR), H-infinity, particle swarm optimization etc [2-6]. 
Some of the examples are described here. In [2] authors have 
shown the control inputs 5 £ and 5 E to provide robust 
performance when compared to the other damping controllers 
by applying a phase compensation control technique with 
respect to state space variable speed. In [3] authors have 
presented iterative particle swarm optimization (IPSO) based 
UPFC controller to achieve improved robust performance and 
to provide superior damping in comparison with the 
conventional particle swarm optimization (CPSO) for the 
control inputs g £ and In [4] author has presented multi 
machine system, where some of the states having larger 
settling time with conventional LQR are well regulated with 
multistage LQR. 

In the current paper, for the modified SMIB linearised 
Phillips-Heffron model, after doing a preliminary control 
analysis with individual inputs and coordinated inputs, a 
switching strategy between individual controllers is 
suggested for UPFC devices such that the steady state 
response and settling time will be moderately better for all 
the four state space variables simultaneously for either cases 
of individual inputs as well as coordinated inputs. 

Paper is organized as follows, in Section II, modified SMIB 
linearised Phillips-Heffron model is described. It is followed 
by some preliminary analysis, first with individual controllers 
later with coordinated controller in Section III. Section IV 
describes the switching model for Philips-Heffron plant with 
UPFC controllers along with the proposed switching rule. 
Results and analysis follow in the concluding section. 

H. DYNAMIC MODEL OF POWER SYSTEM WITH UPFC 

H.F Wang has presented the following state space model 
for the modified SMIB linearised Phillips-Heffron power 
system [1,5]. 

x = Ax(t) + Bu(t> (1) 
Where, the state variables are the rotor angle deviation (Ao), 
speed deviation (A&>"\, q-axis component deviation (AE^), 
field voltage deviation (AEf^') and input variables are 
modulating index and phase angle of shunt inverter 
(m E , fi E "i and modulating index and phase angle of series 
inverter (rr.g. 5g A and B represent the state and control 
input matrices given by 
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All the relevant k-constants and variables along with their 
values used in the experiment are described in the appendix 
section at the end of paper. 

in. PRILEvIINARY OPTIMAL CONTROL ANALYSIS 

In this section, a preliminary analysis is done by 
controlling modulating index of shunt and series inverters 
rQ E and rrig ( me two chosen inputs for the current research) 
using LQR based controllers, in order to gain some insight 
into the system behavior and to arrive at a suitable switching 
strategy. Analysis is done in two stages. In the first stage the 
conventional optimal control(COC) analysis is done by 
selecting rn^ or as tne control inputs individually resulting 
in two separate Single Input Single Output (SISO) systems, 

x = Ax + Bu 
Where g = g the first column of the B matrix for the 

input m E , and B = B 3 for the input m B 
The Control law is given by 

u — r — Kx 

Where, K = K g and K = K j, are the controller gains for 
the inputs ir. E and tn B respectively. Both K s and K b were 
designed by conventional LQR method and state variables 
were analysed. Refer Fig. 1 to 4. In the second stage COC 
analysis is done by selecting both and as the coordinated 
inputs resulting in a Multi Input Multi Output (MIMO) system 
with 

B = [B ± fl 3 ] 

Now controller gain K is 2x4 matrixes for this MIMO model 
obtained by LQR algorithm for MIMO system. Analysis re- 
sults for all the state variables are presented below in Fig. 1 
to 4. 
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Figure 1 : COC for Rotor Angle Deviation 
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Figure 2: COC for Speed Deviation 
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Figure 3: COC for q-axis Component Deviation 
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Figure 4: COC for field voltage Deviation 
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In view of the above optimal control analysis, investigation 
of Fig. 1 to 4 reveals that: 

a) Any of the above COC can not provide better performance 
for all the four state space variables with respect to peak 
overshoot and settling time. 

b) m-£ Provides better performance for q-axis component 
deviation. 

c ) m u Provides better performance for rotor angle and speed 
deviations. 

d) The coordinated inputs and provides better performance 
for the field voltage deviation. 

This optimal control analysis suggests that suitable switching 
between controllers and (corresponding to inputs and) 
may improve the steady state performance of all the four 
state space variables. 

IV. Proposed optimal feed-back switching control 

In this section, mathematical modeling of Phillips-Heffron 
system with UPFC devices as a switched linear systems and 
the proposed switching algorithm will be explained. 

A. Switched Linear System 

Switched systems are composed of a group of 
subsystems guided by a switching law that governs the 
change among the subsystems. Use of appropriate switching 
in control has proved to give better performance when 
compared to the performance of a system without switching 
control. A switched-linear system model (refer Fig. 5) for the 
current problem is as follows: 

x(t) = A m x(f) + Bu (2) 
y(t) = Cx(t) (3) 



Conffioller 1 



CtKBroda : 



Switch 




Plant 


— ► 



Figure 5: General implementation of switched linear systems 

The switching strategy 9(jt\ shown in (2) takes values 1 and 
2 based on switching rule decided by the supervisor leading 

to closed loop A 1 — A — BK ± an d A 2 — A — BK 2 - The 
controller gain vectors can be obtained from, linear quadratic 
regulator theory. For the sake of completeness LQR theory is 
now briefly described [4] . The LQR controller generates the 
parameters of the gain by minimizing the error criteria in (4). 
Consider a linear system characterized by (1) where (A, B) is 
stabilizable. Then the cost index that determines the matrix K 
of the LQR vector is 



/ = -J^(x T Qx + u T Rii)dt 
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Where Q and R are the positive-definite Hermitian or real 
symmetric matrix. From the above equations, 

K = -R~ 1 B I F (5) 
And hence the control law is, 

u(t) = -Kx(t) = -R- 1 B T Px(t) ( 6 ) 
In which P must satisfy the reduced Riccati equation: 

PA+A T P- PBR~ 1 + B T P + Q = 0. (7) 

The LQR function allows you to choose two parameters, R 
and Q, which will balance the relative importance of the input 
and state in the cost function that you are trying to optimize. 
Essentially, the LQR method allows for the control of all 
outputs. 

Here the two controller gains ^ and model the UPFC 

gains with respect to rrij- and nig. The controller gain is the 
primary controller and is the secondary controller. Thus the 
two closed-loop parameters will now be A± = A — BK± 

and A. 2 = A — BK 2- m order for to correspond tom E , 
define K± = [K' g 0]' and for ft^to correspond to mE , 
define K z = [0 K^]'. Now note that, 

Eig(A t ) - EigiA - BK ± ) - EigiA - BK e ) and 
EigiAz) = EigiA- BK Z ) = Eig(A-BK b ~). 

B. Switching Algorithm 

The switching control algorithm based on [7, 8] can be 
explained in the following steps: - 

1 . Define jjjL as the primary controller for ni £ and for m B . 
where A± = A — B K± asymptotically stable and 

A 2 — A — BKi not necessarily stable. 

2. Determine ]~Q.by solving the algebraic Lyapunov Equation 

4r + iv ± = -c T c 

3. Using, A* — A — BK™ define the switching matrix 

s = -(Alr a + r A 2 + c T c) 

4. Now, the switching rule is, use secondary controller ft 2 
with 

0(0 = 2 if < x,sx » 
— 1 otherwise 

Results and conclusions 

The experimental set-up to test the proposed algorithm 
consists of linearised Phillips-Heffron model of SMIB installed 
with UPFC described by A and B (modeling and) matrices 
below. The primary controller and alternate controller are 
obtained by solving Riccati equation using R=l and. The 
matrix C is a vector with zeros along with 1 in anyone position 
depending on the state variables on which the peak overshoot 
and settling time is based. The proposed optimal feed-back 
switching rule S between two controls vector of and is also 
given below. 
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The dynamic response curves for the four state space 
variables rotor angle deviation (A3), speed deviation (AcSji, 
q-axis component deviation (AE{.~), field voltage deviation 
(AE,-£ "i are plotted as shown in the Fig. 6 to 9 with the legend 

Switch Ya-£ an d nig f° r me proposed optimal feed-back 
switching control (OFSC) damping controllers. Inorder to 
show the effectiveness of our proposed method settling time 
is also tabulated for the COC and proposed OFSC. 

From Fig. 1 to 4 of COC and Fig. 6 to 9 of OFSC and Table 
one conclude that the proposed optimal feed-back switching 
control provides robust performance in the steady state 
period and moderately better performance in the settling time 
in all the four state space variables simultaneously compared 
to system response with optimally controlled individual 
inputs without switching as well as optimally controlled 
coordinated input (MIMO model). 




Figure 6: OFSC for the Rotor Angle Deviation 
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Figure 8: OFSC for q-axis component Deviation 
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Figure 9: OFSC for field voltage Deviation 
Table I: Comparison of settling time for COC and OFSC 
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APPENDIX 

Synchronous Machine: 
H = 4.0,U = 0.0, TJ = 5.044 
Excitation System: 
k A = 100,T A = 0.01 

k Constants for the nominal operating conditions: 
fc L = 0.5661, fc : = 0.1712, fc, = 24583 



fc 4 = 04198, fc 5 = -0.513, fc 6 = 0.3516 
fcpa = 0.3795, fe^ = 1.1628, k e „ = -0.4591 
ftpi = 0,1864, k qh = 0.2S55,fc L ,j = -0.1096 
fe pSf = 1.1936, k qSs = -0.0380,^ = -0.031: 
ftpss = 0.0529, = -0.0423, fc^ = 0.0189 



Figure 7: OFSC for Speed Deviation 
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